function up = upsilonp(q,param)

    sigma = param.sigma;
    epsilon = param.epsilon;

    up = (sigma-1)/sigma*q.^(-1/sigma); %(sigma-1)/sigma*exp((1-q.^(epsilon/sigma))/epsilon);

end